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Abstract 

The generalization of Kasteleyn and Fortuin clusters formalism is introduced in XY 
(or more generally 0(n)) models. Clusters geometrical structure may be linked to spin 
physical properties as correlation functions. To investigate percolative characteristics, 
the new cluster definition is analytically explored in one dimension and with Monte Carlo 
simulations in 2D and 3D frustrated and unfrustrated n-clock models. 



1 Introduction 



The idea to describe long range correlations and coherency in spin systems from a geomet- 
rical point of view dates thirty years back. In the late sixties this project was accomplished 
for Ising models by Kasteleyn and Fortuin (KF), who developed a method to give intrinsic 
definitions of clusters of spins which might describe with their percolation characteristics 
the structure of correlations patterns. Cluster formalism and percolation tools [ffl, 
have proven extremely useful in the understanding of critical phenomena of Ising models 
||. [| • Among the many results, very interesting is the discovery of the possibility to de- 
scribe spin correlations trough percolative connectivity functions and the consequent link 
between thermodynamic critical behaviors and cluster fractal structures (see ||). The in- 
dividuation of "physical" clusters of spin introduced within this approach, has been also 
successfully exploited by Swendsen and Wang (SW) || [7| to develop fast Monte Carlo 
(MC) dynamics, for unfrustrated Ising models, based on cluster update and later on to 
drastically improve simulations in frustrated systems too || || . 

Nevertheless, the discussions about the extension of cluster formalism and perco- 
lation concepts to continuous spin systems as XY or 0(n) models is still open, and the 
equivalent of KF cluster in such systems is not known. Wolff [fJJ has proposed a cluster 
definition in XY models based on a smart application of KF rule to spin projections along 
random directions. The clusters so introduced have proven to have a percolative critical 
temperature exactly equal to the thermodynamic one in unfrustrated XY models, also 



if their ultimately connections with spin properties is not understood ||11|| . One of the 
successes of KF approach consists in the clarification of the links between the clusters and 
the physics of the spins. 

In this paper we try to address a possible generalization to XY and 0(n) models 
of KF approach to Ising like systems. This extension leads to new cluster definitions. 
In the spirit of Kasteleyn-Fortuin and Coniglio-Klein (CK) works, we try to focus the 
relations between clusters and spins, disregarding, at the moment, applications to efficient 
MC algorithms. Specifically we try to introduce concepts and tools to manipulate the 
the structure of such clusters in this larger context, as done by KF and CK in Ising 
like systems. In ferromagnetic models, the new clusters generally individuate regions of 
statistically coherent spins, i.e. almost parallel spins, and describe the physics of such 
aggregates. For sake of clarity, before passing to such a generalization (presented in section 
3 and 4), in section 2 a fast outlook to KF original approach to Ising systems is given, 
to define notations and concepts used in what follows. Later, the properties of these 
new clusters are analytically studied in one dimension and via MC simulations in n-clock 
models in 2D and 3D, where the thermodynamic transitions have different properties. 

2 Cluster formalism for Ising spin Hamiltonians 

Let us consider an Ising system of spins Si = ±1 with Hamiltonian: 

m{{Si}) = - E {JuSiSj - \Jiji) (i) 

<i,j> 

where {Si} is the spin configuration, the sum is over all interacting spin pairs and, as 
usually, (3 = 1/ksT. The constants in the Hamiltonian have been chosen for future 
convenience and for simplicity we can consider Jij = J > Vi,j, i.e. we take an Ising 
isotropic ferromagnet. 



The cluster formalism to describe "droplets" of Ising spin || was originally developed 
by Kasteleyn and Fortuin and later on in a different approach by Coniglio and Klein 
@. It is based on the mapping of the original model described by Hamiltonian ([!]) into a 
new model in which spin couplings have infinite or zero strength. The mapping consists 
in stochastically changing the interactions Jjj between spin pairs to new values J[a (to 
define clusters we are interested in the limit \J[j\ — > oo or J[- — ► 0), in such a way that 
the two models, with new and old interactions, are statistically equivalent. 

To this aim, following KF, let's consider then each single couple of interacting spin 
Si and Sj, and suppose J[j G {0, J', —J'} with J' > 0. If the change of to a new 
interaction J[j = J' occurs with a weight po, to a J[j = — J' with p n and to J-j = with 
q, then the sum of statistical weights of a spin configuration {Si, Sj} in the new possible 
models is 0,0]: 

W(S h Sj) = q + p e J '^-V + w-nw+i) (2) 

It is to be noticed that we are slightly modifying KF-CK original approach in which p n 
was a priori set to zero. To impose the statistical equivalence of the original and mapping 
models we must then require that a given spin configuration has the same weight: 

exp(-/3W({^})) = J] W(Si,Sj) (3) 

<i,j> 

So for each pair of interacting spins we must require: 

e J(S tS] -l) = g + pQe J'(S iSj -l) + p ^-j' { S iSj +l) (4) 

where the po, p-n and q are unknown temperature functions. 

To introduce the definition of clusters of spins we must consider the limit J 7 — > oo. 
Two spins Si and Sj connected in the new model by an infinite interaction must have then 
a definite reciprocal direction (i.e. parallel if J[j = +J' and antiparallel if J[j = —J') to 
have a non zero weight, otherwise, if disconnected, they are completely independent. Thus 
in the new model the clusters are naturally defined as the maximal sets of spins connected 
by oo interactions (called bonds). The deletion (J?. = 0) or the freezing (| J[j\ = oo) of the 
original interactions leads to the contraction of the spin lattice in independent fundamental 
units: the clusters. In the limit J' — > oo, eq.(§) becomes: 

^(SiSj-i) = g + po s s . ;S . + p„5 s . ; _ s . (5) 

This is a linear system of two equations with three unknowns q, po and p n , and so it is 
possible to introduce physical constraints to select some definite solution [[J. To this aim 
let's introduce the connectivity function 7^ , which is one if spin Si and Sj belong to the 
same cluster and zero otherwise. It is possible to show that connectivity is always greater 



or equal to spin correlation R 12 



(SiSj) = ( 7 Ji) - {ll) < (4) + (7S> = <7^> (6) 

where 7]^ (7^) is one if Si and Sj belong to the same cluster and are parallel (antiparallel). 
A criterion which proved to be extremely important to select the definitions of interesting 
clusters (i.e. the relative value of q, po and p n ) is to make connectivity as close as possible 
to correlation, i.e. to minimize connectivity as a function of q, po and p n [§, 14]] : 



(jij) minimum 



(7) 



This natural condition, which essentially corresponds to select clusters whose structure 
resembles the correlation patterns in the system, has given excellent results in frustrated 



and unfrustrated Ising spin systems || [LJ]. In the case of the simplest approximation in 
which we consider just a single couple of interacting spin Si and Sj, the mean connectivity 
is (jij) = (p +Pn)/ (1 +e~ 2J ), and to impose condition (^), with constraints < q,Po,PTr, 
naturally leads to p n = (or analogously p = if J < 0). In the present case, this results 
may be also simply obtained by directly imposing (SiSj) = (jij). The solution: 

Ptt = ; p = l-q=l- e~ 2J (8) 

is the well known result by Kasteleyn-Fortuin [[!]] and by Coniglio-Klein ||. Within this 
context it is possible to show that the partition function of a Q-Potts model |4[] may be 
written as (Q = 2 corresponds to Ising model) |], ||: 

Z Q (J)=Y,<l lAl P l o Cl Q N{C) (9) 

a 

where po = 1 — e~® J = 1 — q, N(C) is the number of clusters in the bonds configuration 
C (i.e. the set of oo interactions), \C\ (resp. \A\) is the total number of bonds in C (resp. 
of absent bonds or zero interactions), and J2c * s the sum over a ^ bonds configurations. 
Eq.@ gives the Ising partition function in terms of the partition function of a correlated- 
percolation model [I3|. Moreover, with KF solution, for an Ising ferromagnet, eq.(|5|) 
becomes: 

(SiSj) = (Hj) (10) 



3 Cluster formalism for XY spin Hamiltonians 

Let's examine now the problem of cluster definitions in continuous spin systems. We con- 
sider an XY model, but the same arguments may be extended to 0(n) models. Specifi- 
cally, we consider a system of planar spins, Si, with pair Hamiltonian: 

PHy = -{Jij CO8(0i - dj) - \Jij\) (11) 

where 9i is the phase of spin Si, as above we suppose (3 = l/ZcgT and for clarity = 
J > 0. The constant in the Hamiltonian has been chosen for convenience so that two 
ferromagnetically interacting spin have zero energy when they are parallel. 

Following the idea proposed by KF, we map the original model described by Hamil- 
tonian (|Tl~D into a new model in which the pair Hamiltonian between interacting spins is 
stochastically changed to new functional values H[-, in such a way that the two models 
are statistically equivalent. As above, to individuate clusters we are interested in the limit 
H[- — > oo or H-j — > 0. The main difference with the previous section will consists in the 
fact that many choices for H' are necessary, but the arguments will be the same. Let's so 
define new pair Hamiltonians characterized by a new variable 

(3H{Si, S 3 ; = -(4 cos(0, -6 3 - - |JT.|) + C{J' i3 ) (12) 

where the parameter is - G [0, 2ir], J-j = J' > and C(J') is an adjustable regularization 
function for the limit J' —>■ oo. 



To impose the statistical equivalence of the original and the mapping models, we 
then require that spin configurations have the same weight in both of them, so if we define 
p(4>') as the statistical weight to map the pair Hamiltonian (|TT|) into H(Si, Sj] 0'), and q 
the weight to map it into a zero energy interaction, the equation corresponding to the (Q) 
becomes: 

e-ZHWA) = q + F pU/y-mSi&Mdtf (13) 
Jo 

As above, to define clusters we consider the limit J' — > oo. Two spins Si and Sj 
connected in the new model H(Si, Sj] $ by an infinite interaction, must have a definite 
reciprocal direction (i.e. 0, — = to have a non zero weight, otherwise, if not con- 
nected, they are completely independent. Thus in the new model the clusters are naturally 
defined as the maximal sets of spins connected by oo interactions (bonds). In contrast 
to the Ising case we now have much more than just two kind of bonds (between parallel, 
Po, or antiparallel spins, p w ). Following this method it is then possible to generalize the 
procedure of deletion (J[j = 0) and freezing (J^ = oo) of the original interactions. Also 
in this case different clusters are independent (if q is a function of 0^ — 9j or if the infinite 
limit is not systematically taken, then one has interacting clusters). In the infinite limit 
J' — > oo, eq . ([13|) becomes: 

e-PH^SuSj) = q+ f 27T p^sy. _ _ = q + p (0. _ ) (14) 

Jo 

where the function C(J') has been absorbed to regularize the definition of the 5-function 
in the interval [0, 2tt] with argument defined modulus 2n. Eqs.flU]) is a linear functional 
equations in the unknown functions p(4>', Jij) and q(Jij). 

A solution of eq.(|l4]) is suggested by the reasonable limit behavior p(<f)') — > if 
J — > 0, or alternatively by the condition of local minimal connectivity (see below): 

q = e - 2J ; p{(j)') = eJ^W-i) - q (15) 

This solution reproduces in the Ising case the results by Kasteleyn and Fortuin given in 
eq.®. 

The clusters are operatively individuated by the conditioned probabilities: 

p^'\0. _ 0.) = - b, - 0V H(SiA) ; q{9i - 0j) = q ■ e? H ^ . (ie) 

which may be respectively interpreted as the conditioned probability to substitute the 
original interaction in the mapping model with a bond of the kind <p' or with a zero inter- 
action, given the spin configuration {Si, Sj} (note that this probabilities are completely 
independent on the choice of the constants for the energy of the ground state). These 
conditioned probabilities may be used to implement MC cluster algorithms because they 
contain the necessary information to build clusters from spin configurations and it may 
be proved that algorithms based on these probabilities satisfy detailed balance principle 
II [J. They are the generalization to XY of KF bond conditioned probabilities in Ising 
systems. 

Eq.fllij) may also be considered directly as the starting point to define clusters of 
bonds variables {0^}, avoiding at all to introduce the procedures of Hamiltonian mapping 
and definitions (0). In this perspective eq. flHD is just a way to introduce a statistical 
systems of variables of spin and bonds ({0i}, {4>ij}) with the following peculiar proper- 



ties 151: the marginal distribution of the {OA is exactly equal to the Boltzman weight 



e -/?ff({0i})| the conditional distribution of the {0^}, given the {Oi}, is exactly expressed 
by eqs.(|16|); the conditional distribution of the {Oi}, given the {0^}, correspond to the 
above given definition of clusters, i.e. two interacting spin Si and Sj belonging to the 
same cluster must have the definite reciprocal direction 0, — 9j = 0^-, otherwise, if discon- 
nected, they are completely independent. Note that only the sets {0^} such that given 
any two sites h and k the quantity A9 = YJh 4>'ij is independent of the "integration" path, 
are allowed. 



4 Relations between thermodynamics and percola- 
tion 

The previous section was devoted to introduce a simple generalization of Kasteleyn- 
Fortuin and Coniglio-Klein clusters in XY models. Now we face the problem to work 
out some main relations between percolative and thermodynamic quantities. Easy exten- 
sions may be given for general 0(n) models. 

The partition function of the XY model, from eq. fllil) , may be written as: 

Z = II r dOie-^<^> H{S ^ ] = j2q lAl P(C)(27i) N ^ (17) 

* c 

where \A\ is the number of absent bonds on the lattice fixed the bonds configuration C, 
N(C) is the total number of clusters in C, and by definition P(C) = I\ < ij >e cP(4 ) ij)- The 
sum, J2*c> is intended over all possible bonds configurations, C, (note that two bonds 
configurations are distinguished by their geometry and by the kind of bonds {0^} they 
have), and specifically: 

e=e/( n (is) 

C* C <ij>£C 

with n<jj> 6 c the product over all bonds present in the configuration C, {0^} the set of 
indexes of such present bonds and S(C, {4>'ij}) a function nonzero only if the configuration 
C and the set are compatible (i.e. if the sum of 0^- between whatever fixed extrema 
h and k along a chain of present bonds, is independent of the path, i.e. clusters are well 
defined because such quantity is exactly the phase difference between Sh and Sk)- 

The percolative quantity to be compared to the thermodynamic two point correla- 
tion function is the pair connectivity c(i,j), defined as: 

c(hj) = / c(i,j,<p)d<p (19) 

J 

where c(i,j, 0) = (jij(4>)) and 7r/(0) = jij ■ S(9i — 9j — 0) which is zero if spin Si and Sj 
do not belong to the same cluster or have a phase difference 9, L — 9j ^ 0. c(i,j,<f>) is so 
the probability of spin i and j to belong to the same cluster with a phase difference 0. 
It is possible to show that the pair correlation function g(i,j) = (Si ■ Sj) is given 

by: 

g(i,j)= cos(0)c(i, j,0)d0 (20) 

.In 



Eqs.([T9|) and ( pO|) imply g(i,j) < c(i,j), analogously to Ising systems where eq. (||) holds. 
A consequence of this proven inequality is that T c < T p , where T c and T p are defined as 
the temperatures where respectively the magnetic susceptivity x an d the mean cluster 
size S = J2 s n s s2 EH ( n s is the number of clusters of size s and J2 S is the sum over just 
finite clusters) become singular. These relations naturally suggest, in analogy with the 
Ising spin case, the criterion of minimal connectivity to select clusters definitions, i.e. to 
impose the condition: 

c(i,j) — > minimum (21) 

where c(i,j) has to be minimized respect to p(<p). As anticipated above, imposing eq. (|21~D 
just for each single couple of interacting spin Si and Sj, directly leads to select the 
solution given in eq.(|l^) for eq.([L4]). It is to be noticed that this solution is however 
just the simplest extensions of KF result. In facts more general solutions must be found 
as showed below, but the general tools introduced to link percolative and spin properties, 
allows to exploit, for XY models, the many techniques to individuate and manipulate 
clusters known in the literature for Ising systems (see |], §, |], 0])- 111 what follows we 
will restrict however to consider the simple solution given in eq . (|l~5|) . 

The relations above reported indicate that thermodynamic spin quantities may be 
generally expressed in terms of cluster properties. For example it is possible to link the 
mean energy E with geometrical quantities. In the case of an isotropic ferromagnetic XY 
model it results: 

E/N = -<7oifyl*(p(0oi)/<7)> + ^ln(g) (22) 

where and 1 are two of the N interacting pairs of spin in the system. As eq. (|l7"D 
is the natural generalization of eq.(|9|), so eq.(^) is the extension of the corresponding 
energy-bond relation in Ising systems [0, |T6|j. 



5 Clusters in one dimensional XY model 

To understand the properties of the above defined clusters it may be interesting to an- 
alyze the question in some details. In a one dimensional XY model of nearest neighbor 
interacting spin |T7||, the geometry of the above defined clusters corresponds to chains of 
bonds, and the problem is extremely simplified. So it is possible to prove that, adopting 
solution (pT5|), the partition function, eq. (|17D , in the case of an isotropic XY ferromagnetic 
chain, is: 

Z 1D (J)=J2l lAl P lCl W NiC) (23) 

c 

where J2c ^ s J us t the sum over all graphs of bonds on the chain, P = 2ire~ J (Iq(J) — e~ J ) 
(I n (x) is the imaginary argument Bessel function of order n). From eq. (|23|) it is possible 
to see that the partition function of an XY chain may be written as that of a Q-Potts 
linear model, Zq (see eq.(|9])), times a simple factor. Specifically: 

Z 1D (J) = K ■ Z Q (J Q ) (24) 

where N is the total number of interactions, Q = 2ir, Q ■ Jq(J) = hi(l + P/q) (Jq ~ J if 
J — > oo or J — > 0), and K(N, J) = e - N ( 2J -Q J Q) _ 

We are concerned with clusters and spin properties, and in the context of the linear 
model it is possible to prove a definite relation between correlation g and connectivity c: 



g Uj) = e r/i c(i,j) 



(25) 



where r = \i — j\ is the number of spin between Si and Sj plus one, and 

C 1 (J) = H(W)-e~ J )/h(J)} (26) 

Moreover the connectivity c(i,j) at temperature T = 1/J for the XY model, is equal 
to the connectivity CQ(i,j) of the KF (or CK) clusters introduced above in the Q-Potts 
model with T = l/Jg(J) and Q = 2n: 

c(i,j)\j = c Q (i,j)\j Q{J ) (27) 

These results should shed some light on the connections between cluster connectivity and 
spin correlation in one dimensional XY models. A trivial consequence of all these relations 
in ID is that T v = T c = (it would be hard to find clusters with a bit of randomness with 
T p > in one dimension), but this coincidence does not hold for the critical behavior. 
Defining £xr as the correlation length in the XY model and ^Q-p tts a s the mean cluster 
radius in its Q-Potts equivalent, eq. (p5|) imposes ^ x \- = £,nlp otts + but, in ID, at low 



temperature ^xy ~ T 1 while ^Q~p tts ~ e A ^l T JI3L [T7|], and so clusters quite loosely 
express spin-spin correlations. 



6 MC results for XY model in higher dimensions 

The analytical problem concerning the structure of clusters in XY models in higher 
dimensions, is, worthless to say, much more difficult. We present then some Monte Carlo 
results about clusters properties (defined from solution (|15|)), in two and three dimensions. 

MC simulations were done using a standard Metropolis spin flip algorithm [jXS | on 
n-clock models on a square or cubic lattice described by Hamiltonian (|TT|) (in the following 
if not specified we will consider the isotropic case Jjj = J > 0), whose spin Si have a 
phase Qi = 2nm/n with m £ {0, n — 1}. 

Let's briefly examine our MC results in two dimension. In the case n = 2 we exactly 
recover the 2D ferromagnetic Ising model. Our MC simulation indicate the well known 
result of equal critical temperatures T c = T p = 2.269 (all temperature are measured in unit 
of J) and a percolative critical phenomena characterized by Ising exponents. Moreover, 
the MC dynamic based on the above defined clusters is just the Swendsen-Wang dynamic 
and the phenomena of critical slowing down drastically reduced || 0. 

For n > 4 such correspondence is no longer verified. Such a result may be expected 
because in these cases connectivity and correlation are not coinciding as shown by the 
simple example of just two interacting spin with n = 5 (the cases n = 3,4 may be 
successfully faced with some tricks, in resemblance of the possibility to map n = 3, 4 
clock- models in an equivalent Q-Potts [|J). 

The percolation critical temperature decreases for increasing n, and approaches a 
plateau in the large n limit. In actual facts, via MC simulations, for n = 36, we find that 
the percolation point is at 

T p = 1.69 ±0.03 



to be compared with the 2D XY critical temperature at T c ~ 0.89 |7|, |TT|. We find that 
the percolation critical exponents are in the universality class of random percolation as 
expected because there is no thermodynamic transition underlying the percolative one: 
the critical exponents, measured via a finite size scaling analysis fl8j (reported in Fig. 



1), are u = 1.33 ± 0.05 and 7/1/ = 1.79 ± 0.05 in perfect agreement with 2D random 
percolation exact values v = 4/3 and j/u = 43/24. These exponents, in percolation, 
characterize respectively the divergence of mean cluster radius £ and mean cluster size S: 
and 5* 



IT — T n 



\T-T n 



/' " ' * i- 1 -"-pi 0- 

This behavior is observed, as may be easily suspected, in frustrated or disordered 

systems too. We tested, via MC, the Fully Frustrated XY model (FF) and the ±J XY 

Spin Glass (SG) (see references in ]19[), where we found the same percolative critical 

exponents and (see the scaling analysis in Fig. 2 and 3): 



T, 



FF 



1.61 ± 0.03 and T/ G 



1.64 ±0.03 



This value for the percolation transition in the FF model is above the critical region 



located around T ~ 0.4 0.5 (see [19]). Also Wolff's clusters show a percolation point 
well above the critical region |TT|, but it is possible to introduce their direct generalizations 
whose T p may be pushed closer and closer to it |2tJ . 



The same kind of results are found in three dimension. For n = 2 we recover the 
well known properties of KF or SW clusters in the 3D Ising model T p = T c ~ 4.5 and 
v ~ 0.62, j/u ~ 1.97 (see Our MC runs show, for n = 36, T p = 3.75 ± 0.05 and 

3D random percolation critical behavior with v = 0.87 ± .05 and 7/1/ = 2.00 ± .05 (see 
Fig. 4). These values are to be compared to the results of 3D XY T c ~ 2.2 and v ~ 0.66, 
7/1/ ~ 1.98 (see [2"2|). Essentially the same values are found for the 3D ±J XY Spin 
Glass. 

It is interesting to note that, as expected and discussed above (see also [§, |i~4"|), 
whenever the gap between T p and T c becomes finite, SW like cluster algorithms for MC 
simulations become unable to reduce critical slowing down. 

After these MC results, the panorama we get illustrates that the straight generaliza- 
tion to XY models of KF clusters, given in eq. (p~5|) , has not the peculiar properties of KF 
clusters in Ising like systems: the thermodynamic and percolative transition are no longer 
coincident. In Ising models more complex procedures have been introduced to individuate 
physical clusters, as those proposed in HP, It would be interesting to verify if the 
extensions to XY of such procedures according the lines proposed above, have the same 
percolation properties here found or new interesting results can be obtained. 



7 Summary and conclusions 

In analogy to Kasteleyn and Fortuin and Coniglio and Klein works, cluster of nearest 
neighbor spin in XY models may be defined as the sets of spin connected by bonds 
according definite rules. The clusters divide the original lattice into independent regions 
of statistically coherent spins. Kasteleyn and Fortuin percolative concepts and tools to 
link clusters and spin properties, which proved to be so useful in Ising systems, can be 
so extended to XY models. In these models, at a first simple level, KF clusters may be 
defined by separately looking at just each single couple of interacting spins. Consequently, 
bonds are introduced between them according the definite probability distribution given 
in equation (|i~5|) (as a matter of fact, this is for XY models, but the analog for general 
0(n) is absolutely similar). 

Nevertheless, at the simple level here explored, many differences appear with Ising 
systems. It is known that in the Ising model these clusters have a percolation point 
which, imposing condition (|2~1"D, may be pushed to coincide with the critical one T c , and 



has percolative exponents in the Ising universality class ][l fj. The properties of the new 
clusters may be studied analytically for continuous spin systems in one dimension. Here 
eq.(p5|) implies that clusters mean square radius, differently from Ising KF case, is no 
longer coincident with the correlation length in the system. Numerical results show the 
same behavior in higher dimensionality, where moreover cluster percolation point T p is 
different from T c . Phenomenologically, the temperature T p is the point where regions of 
almost parallel spin (in ferromagnetic models), i.e. regions of coherent spins, percolate in 
the system. The transition corresponding to this point, for clusters defined by eq.fllSf)) is 
in the random percolation universality class. 

The occurrence of a finite gap between T p and T c is found in Ising spin systems 
when frustration is present. In this general criterion to close such a gap has been 

proposed || ||, [TJj]. Exploiting the results here presented, it is possible to apply such 
a criterion to frustrated and unfrustrated 0(n) models too, and in perspective give a 
percolative description of their critical behaviors in analogy to the known results for Ising 
like models. This approach would lead to a change of the bond probability distribution 
given in eq. fllBl) . In unfrustrated and frustrated Ising models a definite physical origin 
has been associated to the percolation point T p |[23j| . It is then natural to speculate on it 
in continuous spin models too. 

The criterion introduced in || |], [14j] is actually suited to develop efficient MC cluster 
algorithms in Ising systems. The perspective to go further in such a direction also for 
0{n) models, is very appealing. 

The author is grateful to Prof. Antonio Coniglio for stimulating discussions and 
suggestions. 
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Figure Captions. 

Fig. 1. Finite size scaling of the mean cluster size S of clusters defined by bond prob- 
abilities given in eq. fllB"!) , in the 2D XY ferromagnet. The scaling parameters are 
T p = 1.69 ± .03 and v = 1.33 ± .05 7/1/ = 1.79 ± .05 

Fig. 2. Finite size scaling of the mean cluster size S of clusters defined by bond prob- 
abilities given in eq.(|16|), in the 2D XY Fully Frustrated. The scaling parameters 
are T p = 1.61 ± .03 and v = 1.33 ± .05 j/v = 1.79 ± .05 

Fig. 3. Finite size scaling of the mean cluster size S of clusters defined by bond proba- 
bilities given in eq.(|l6|), in the 2D XY ±J Spin Glass. The scaling parameters are 
T p = 1.64 ± .03 and v = 1.33 ± .05 j/u = 1.79 ± .05 



4. Finite size scaling of the mean cluster size S of clusters defined by bond prob- 
abilities given in eq. (p~6[) , in the 3D XY ferromagnet. The scaling parameters are 
T p = 3.75 ± .05 and v = 0.87 ± .05 7/z/ = 2.00 ± .05 
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